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ABSTRACT 


A  model  has  been  developed  to  calculate  the  steady  state  vibrational 
distribution  of  an  anharmonic  oscillator  under  conditions  of  thermal  non- 
equilibrium.  The  model  includes  vibration-vibration  and  vibration- 
translation  collisional  processes  as  well  as  radiative  decay.  Numerical 
calculations  are  shown  to  be  in  reasonably  good  agreement  with  experi¬ 
mental  data  obtained  in  electrically  excited  N^-CO  mixing  systems.  It 
is  shown  that  radiative  decay  may  be  the  dominant  loss  mechanism  in  low 
pressure  ^-CO  experiments. 
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I.  INTRODUCTION 


The  process  of  vibrational  relaxation  under  nonequilibrium 
conditions,  such  that  the  vibration  mode  is  out  of  equilibrium  with  the 
other  degrees  of  freedom,  is  of  considerable  importance  in  many  physi¬ 
cal  phenomena.  The  upper  atmosphere  provides  an  example  wherein 
the  nitrogen  vibrational  energy  is  thought  to  be  determined  by  collisions 
with  free  electrons  and  metastable  oxygen  atoms  and  is  out  of  equili¬ 
brium  with  the  low  translational  temperature.  A  similar  situation 
arises  in  expanding  gas  flows,  molecular  lasers,  and  many  chemical 
reactions  where  the  vibrational  temperature  may  be  much  higher  than 
the  translational  temperature.  In  the  past,  analyses  of  such  nonequili¬ 
brium  conditions  have  been  characterized  by  a  Boltzmann  distribution 
corresponding  to  a  specific  temperature  for  each  degree  of  freedom. 

Recently,  several  theoretical  analyses  have  been  developed  to  avoid  the 
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restriction  of  a  Boltzmann  distribution  in  the  vibrational  mode.  ’ 

There  are  two  basic  colli sional  phenomena  in  the  relaxation  of  a 
diatomic  molecular  gas  from  a  state  of  vibrational  nonequilibrium.  These 
are  the  processes  of  vibration- vibration,  V-V,  and  vibration  to  trans¬ 
lation,  V-T,  energy  exchange.  In  the  following  analysis  we  will  assume 
that  the  translational  and  rotational  degrees  of  freedom  are  in  equilibrium 
so  that  vibration- rotation,  V-R,  or  vibration- rotation-translation,  V-R-T, 
energy  exchange  processes  are  indistinguishable  from  V-T  processes. 

In  general,  the  V-V  relaxation  rates  are  much  faster  than  the  V-T  rates, 


and  under  suitable  conditions  a  gas  may  relax  to  a  state  of  vibrational 
quasi  equilibrium  before  V-T  processes  become  important.  The  colli- 
sional  exchanges  may  occur  between  different  species,  intermode  ex¬ 
change,  or  between  different  vibrational  levels  of  the  same  species, 
intramode  exchange.  For  an  anharmonic  oscillator  the  intramode  ex¬ 
change  may  be  the  controlling  mechanism  in  the  relaxation  process. 

1 

Treanor  et  al  have  considered  the  relaxation  behavior  of  anhar¬ 
monic  oscillators  when  the  vibrational  mode  is  out  of  equilibrium  with 
the  translational  mode  under  conditions  of  V-V  equilibrium  with  no  quanta 
loss  mechanism.  For  these  conditions  they  deduced  non-Boltzmann 
distributions  in  the  anharmonic  oscillator  with  the  vibrational  energy 
pumped  to  the  closely  spaced  upper  vibrational  levels.  These  theoretical 
distributions  are  unrealistic  in  the  higher  levels  in  that  they  neglect 
energy  loss  processes  which  prevent  large  accumulation  of  energy  in 
these  levels.  Bray^has  treated  the  relaxation  process  allowing  for  V-T 
energy  exchange  and  has  obtained  numerical  results  for  quasi  steady 
state  distributions  in  a  single  species.  In  these  calculations  the  loss  of 
vibrational  energy  via  V-T  processes  prevents  the  large  total  inversions 
predicted  by  Treanor  et  al.  These  analyses  of  Treanor  and  Bray  did 
not  include  the  effects  of  radiative  decay  which  can  be  important  in  low 
pressure  systems  where  the  effective  time  between  energy  exchange 
collisions  becomes  comparable  with  the  radiative  lifetime.  Some  pre¬ 
liminary  calculations  have  shown  that  the  radiative  decay  may  be  the 
dominant  loss  mechanism  in  the  ^-CO  mixing  laser.  ^ 
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The  present  paper  describes  steady  state  solutions  for  conditions 
of  vibrational  nonequilibrium  in  gas  systems  composed  of  diatomic 
molecules.  Both  radiative  der.ay  and  V-T  collisional  exchange  are  in¬ 
cluded  in  the  analysis  which  is  applied  to  pure  CO  and  ^-CO  mixtures. 
Vibrational  distribution  functions  are  derived  for  a  variety  of  conditions, 
and  comparisons  are  made  with  existing  experimental  data  in  the  elec¬ 
trically  excited  ^-CO  system. 
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! 
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II.  KINETICS 


Consider  a  mixture  of  two  diatomic  molecular  gases  N  and  M. 
We  assume  that  only  single  quantum  transitions  occur  in  the  collisional 
exchange  processes  of  vibration- vibration  (V-V)  and  vibration-tran- 
lation  (V-T)  as  well  as  in  the  radiative  decay.  The  intramolecular, 
intermolecular  and  radiative  processes  for  species  N  are 


N  +  N  .  ~  N  +  ,  +  N  ,  _  . 

V  v'  vl  1  V1  4.  1 

intra.  V-V 

(1) 

N  +  M  „  32  N  .  .  +  M  , 

v  v"  v+  1  v"4  1 

inte  r.  V  -  V 

(&) 

NV  +  Nv-  ~  Nv+  1  +  Nv 

intra.  V-T 

(3) 

Nv  +  Mv„  *  Nv+  j  +  Mv„ 

inte  r.  V-  T 

(4) 

N  —  N  ,  +hv 

V  v- 1 

rad. 

(5) 

where  the  vibrational  level  is  denoted  by  subscripts  v,  v',  and  v".  A 
similar  set  of  reactions  applies  to  species  M.  One  can  write  the 
master  vibrational  rate  equation  for  the  density  of  molecule  N  in  the 
vth  vibrational  level 
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dN 

v 

dt 
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1  V-T'  v  1  V-T' 
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N  , 
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+  Z 


N-M 


v-  l,v 
P  N 
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v+l,v  /  v,  v-1  v,v+l\  j 

,  +  P  N  -  (P  +  P  N 

1  V-T"  V+1  \  V-T"  V-T"/  v\ 


+  N  ,  /  -  N  /r 
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Z  M  „ 


(6) 


where  Z^  ^  and  Z^_j^  are  the  gas  kinetic  collision  frequencies  of  species 

N  colliding  with  species  N  and  M  respectively,  at  unit  concentrations. 

v.vt  1  v,vtl 

The  probabilities  P  and  P  -refer  to  the  V-V  and  V-T  exchange 

v',v'+l  V-T' 

respectively  and  r  is  the  radiative  lifetime  of  level  v.  The  forward  and 
v 

reverse  rate  constants  are  related  through  detailed  balancing  and  for  a 
Maxwellian  velocity  distribution  corresponding  to  the  translational 


temperature  T  we  obtain; 

v,v+l  v+l,v 

P  =  P  exp, 

v'.v'-l  v'-l,v' 


(K 


Ev'.l'+'VEV+1» 


)]  At| 


(7) 
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th 

E  is  the  energy  of  the  vl  vibrational  level  which  is  approximately 


E  = 
v 


'[ 


(v  +  j)  w  -(v  +  i)2u 


e  e 


] 


he 


(8) 


and  the  molecular  constants  and  wgxe  are  given  in  Ref.  4. 

The  time  rate  of  change  of  vibrational  energy  can  be  obtained  by 
multiplying  the  "master  equation"  (6)  by  v  and  summing  over  all  values 
of  v.  In  the  case  of  harmonic  oscillators,  the  energy  relaxation  equation 
is  considerably  simplified.  For  a  mixture  of  gases  N  and  M  it  can  be 
shown  that 


TV-T 


0m/t 


(9) 


kTV-VGM 


em  <en-en>  "-'8n/t>  -  en,em-em,(1-''Gn/t> 


with  a  similar  equation  for  species  M  where  use  has  been  made  of  the 

5 

Landau-Teller  selection  rules.  Here  is  the  mole  fraction  of  species 
M,  E  is  the  vibrational  energy  per  molecule,  E  is  the  equilibrium  vibra¬ 
tional  energy  at  the  translational  temperature  T,  and  9  is  the  character¬ 
istic  vibrational  temperature,  i.  e.  (Ej-Eg)/k.  Ty  ,p  and  Ty  y  are  the 
vibrational  relaxation  times  for  V-T  and  V-V  processes.  For  times  t 
such  that  Ty_y«  t  «  Ty_rp,  7"ra{j  a  quasi  equilibrium  of  the  V-V  process 
is  obtained.  The  bracketed  expression  in  the  above  equation  is  then  zero 
so  that 


eN 

T 

VN 


0M 

T 

1  VM 


%  - 
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(10) 
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where  Ty*^  and  are  the  vibrational  temperatures  of  species  N  and  M 

respectively.  It  is  clear  that  these  temperatures  must  be  different  under 
conditions  of  thermal  nonequilibrium  and  for  T<  Ty^,  Tyj^  the  species 
with  the  lower  characteristic  temperature  will  have  the  highei  vibrational 
temperature.  This  has  been  previously  derived  by  Teare^  and  indepen-, 
dently  by  Treanor.^  In  the  latter  analysis,  Treanor  obtained  an  analytic 
solution  to  the  master  equations  for  anharmonic  oscillators  allowing  only 
V-V  exchanges.  Their  analysis  exhibits  a  similar  pumping  phenomenon 
between  different  levels  in  the  oscillator. 

In  the  present  paper,  steady  state  solutions  to  the  master  equations 
are  obtained  numerically  for  conditions  of  thermal  nonequilibrium.  The 
equations  considered  are  restricted  to  a  finite  number  of  vibrational  levels 
whose  energy  cannot  exceed  the  dissociation  limit.  The  vibrational  distribu¬ 
tion  in  species  N  is  determined  by  supplying  a  fixed  source  which  commun¬ 
icates  with  all  vibrational  levels.  Under  many  situations  it  is  sufficient 
to  specify  a  fixed  temperature  between  the  first  two  vibrational  levels  of 
species  N.  This  is  equivalent  to  the  inclusion  of  a  source  term  which  only 
communicates  with  the  first  and  zeroth  levels.  Under  these  conditions, 
the  master  equations  may  be  written  in  the  form 


dNo 

dt 


-F0,  1  + 


D,  -  K 


(ID 


dNj 

dt  =  F0,  1  "  Fl,2 


Dj  +  D2  +  K 
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where  K  is  the  net  rate  of  quanta  supplied  by  the  source.  In  these  equations 

F  ,  ,  is  the  particle  flux  between  levels  v  and  v+1 
v,v+l  r 


F 


v,v+l 


(v.v+l 

P  N  N.  - 
...  vi 
1, 1-1 


v+l,v 

p 

i-l.i 


N 


v+1 


(12) 


and  is  the  net  quanta  loss  between  levels  v+1  and  v 

[v+l,v  v,v+l  ] 

PV.T  Nv+1  -  PV.T  NvJ  +  Nv+A+1 


(13) 


These  equations  in  combination  with  the  remaining  master  equations  yield 
the  steady  state  identities 


v>l 


Fv,  v+1  =  Dv+1 


(14) 


and  it  is  clear  that 

K  *  E  Dv+1  (15) 

v=0 


We  assume  a  reflecting  boundary  so  that  there  is  no  transport  of  quanta 
beyond  the  last  level. 

The  numerical  solutions  to  the  master  equations  depend  upon  the 
choice  of  rate  constants.  This  is  in  contrast  to  the  case  of  pure  V-V 
equilibrium  treated  by  Treanor  et  al*  in  which  the  vibrational  distribution 
is  determined  solely  by  detailed  balancing.  Unfortunately,  there  are  few 
experimental  data  on  the  rate  constants  for  V-V  and  V-T  collisional  ex¬ 
change,  particularly  in  the  high  vibrational  levels,  and  one  must  rely  on 
theoretical  calculations.  Since  we  are  interested  in  comparisons  of  the 
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distribution  function  calculations  with  measurements  in  ^-CO  mixtures, 
we  have  restricted  ourselves  to  collisional  exchange  probabilities  for  those 
gases.  The  probability  distributions  used  for  CO  are  shown  in  Fig.  1. 

The  V-T  probabilities  were  derived  from  experimental  measurements  by 

g 

Millikan  and  White  for  the  l-*-0  transition  and  were  applied  to  all  levels 

9 

using  a  correction  for  anharmonicity  based  on  the  S.  S.  H.  model. 

This  model  predicts  that  the  probability  varies  approximately  as 
2/3 

exp  (const.  AE  '  )  where  AE  is  the  energy  spacing  between  the  adjacent 

vibrational  levels. 

For  the  probability  of  intramolecular  V-V  exchange  in  CO,  a 
combination  of  two  theoretical  models  was  used.  The  first,  developed 
by  Sharma  and  Brau,  ^used  a  long  range  attractive  dipole -dipole  potential 
which  is  important  in  the  exchange  probability  in  near  resonant  collisions. 
The  second,  due  to  Rapp,  ^employed  a  short  range  repulsive  potential 
for  the  nonresonant  collisions. 

Because  of  the  large  nonresonance  in  the  first  levels  of  ^  and  CO, 

188  cm" only  the  short  range  calculation  was  used  with  an  assumed 

steric  factor  of  l/3.  During  the  course  of  this  work,  some  experimental 

data  for  the  V-V  deactivation  of  CO  in  vibrational  levels  4  to  1  3  were 

12 

obtained  by  Hancock  and  Smith  in  a  fluorescence  experiment.  The 

intramolecular  data  for  CO  agree  to  better  than  a  factor  of  2  with  our 

theoretical  estimates  as  shown  in  Fig.  1.  Hancock  and  Smith  also 

13 

observed  deactivation  of  CO  by  N2,  and  an  extrapolation  of  their  experi¬ 
mental  probabilities  would  be  less  than  a  factor  of  two  lower  than  our 
assumed  value.  This  is  in  reasonable  agreement  with  other  recent 
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Fig.  1  Collisional  exchange  probabilities  for  both  V-V  and  V-T 

processes  in  CO  at  two  translational  temperatures.  The  data 
points  are  experimental  measurements  of  Ref.  12.  Also  shown 
is  the  radiative  transition  probability  for  CO. 
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RADIATIVE  TRANSITION  PROBABILITY  (secs"1  ) 


14  15 

measurements,  *  Also  shown  in  Fig.  1  is  the  radiative  transition 

probability  for  CO,  These  values  are  based  on  the  measured  v=l  to  0 

probability  and  the  matrix  element  calculation  of  Young  and  Eachus.^ 

_  3 

The  variation  a  v  (Ev~Ev  is  a  close  approximation  to  their 

calculations  for  both  a  linear  and  a  cubic  dipole  moment. 
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III.  RESULTS  AND  DISCUSSION 


In  this  section  we  present  detailed  numerical  calculations  for 
steady  state  distributions  over  a  wide  range  of  energy  deactivation  con¬ 
ditions.  Results  are  first  presented  for  vibrational  distribution  in  pure 
CO  and  then  for  N^-CO  mixtures,  and  finally  comparison  is  made  with 
experimental  data  in  N^-CO  experiments.  In  most  cases,  the  vibrational 
temperature  of  the  lowest  levels,  which  we  designate  by 
9^=  (E^-EQ)/(kf  n  Nq/Nj  )  was  chosen  to  be  2700°K.  This  value  is 
typical  of  the  experimental  conditions. 

A.  V-T  Deactivation  in  Pure  CO 

The  fundamental  characteristics  of  distribution  functions  derived 
in  the  absence  of  radiative  decay  are  of  interest  not  only  in  high  pressure 
CO,  where  the  effects  of  radiative  loss  are  negligible,  but  also  are 
representative  of  solutions  for  homonuclear  molecules  in  general. 

Several  vibrational  distributions  are  shown  in  Fig.  2  for  a  value 
of  9j  of  2700°K.  The  different  curves  were  calculated  for  the  V-V 
probability  distribution  described  in  the  previous  section  but  with  a 

+  2 

variation  of  10-  in  the  magnitude  of  the  V-T  probability  distribution. 
Also  shown  for  comparison  is  the  Treanor  solution  to  the  master  equation 
for  V-V  equilibrium  without  deactivation.  This  latter  result  dependes 
solely  on  detailed  balancing  and  is  described  by* 


N, 


N 


exp 


v+1 


/ErEo  ,  hc\ 
\  kO*  6  e  kT  J 


(16) 
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VIBRATIONAL  POPULATION  (NORMALIZED) 


Tig.  2  CO  vibrational  distribution  functions,  neglecting  radiative  decay, 
for  Qj*  =  2700°K  and  T  =  325°K  for  various  values  of  the  parame¬ 
ter  r,  the  ratio  of  the  V-T  probability  used  in  the  calculation 
to  that  shown  in  Fig.  1.  The  Treanor  distribution  corresponds 
to  r  =  0. 
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V-T  collisions  dominate  V-V  exchange  in  the  upper  vibrational  levels 
because  of  the  large  V-T  probabilities  in  these  levels.  This  upper  level 
deactivation  prevents  the  storage  of  vibrational  energy  and  precludes 
the  artificially  high  inversion  predicted  by  Eq.  (16)  in  high  vibrational 
le  vels. 

At  low  vibrational  levels  where  the  V-T  probabilities  are  insignifi¬ 
cant  relative  to  the  V-V  probabilities,  the  distribution  follows  that  of 
Treanor  et  al  up  to  the  minimum  of  the  latter  distribution.  Beyond  this 
point,  it  can  be  seen  in  Fig.  2  that  the  distribution  behaves  as  vN^  = 
constant  up  to  vibrational  levels  where  V-T  deactivation  is  important. 

From  that  point  on,  the  distribution  reverts  to  a  Boltzmann  distribution 
at  the  translational  temperature. 

The  characteristic  shape  of  the  distribution  function  can  be 
explained  quantitatively.  The  constant  vNv  behavior  is  due  to  the  fact 
that,  for  the  conditions  and  V-V  probabilities  chosen,  near  resonant 
collisions  are  far  more  effective  in  the  transfer  of  vibrational  quanta 
than  nonre sonant  collisions.  Thus,  energy  must  be  transferred  up  the 
vibrational  ladder  by  near  resonant  collisions  rather  than  go  directly  from 
the  source  in  the  first  vibrational  level  to  the  deactivation  region.  Beyond 
the  region  of  the  Treanor  minimum,  the  quanta  flux  supplied  by  the  source 
can  only  be  delivered  from  level  to  level  by  collisions  with  molecules  in 
lower  levels,  i.e.  for  any  level  v  +  1  beyond  the  Treanor  minimum: 

v 

£  Pv,v+1  NvN.«  K  (17) 

i,  i- 1 

i=  1 

where  the  dominant  contribution  to  the  sum  comes  from  levels  close  to  v. 
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By  the  use  of  the  Landau-Teller  relations  it  follows  that  this  is  satisfied 
by 


N  =  constant 
v  \ 


0,T>) 


(18) 


in  agreement  with  the  numerical  calculations.  This  result  is  only  valid 

up  to  the  deactivation  region.  The  region  of  strong  V-T  deactivation 

occurs  when  the  net  particle  flux  cannot  maintain  the  constant  vN^  slope 

to  supply  the  local  deactivation.  A  representative  term  in  the  flux  then 
v,  v+ 1  2  v+ 1» v 

becomes  P  .  N  -  P  ,  N  ,  ,  N  ,  and  by  use  of  the  Landau- 
v,v- 1  v  v-l,v  v+1  v- 1  7 

Teller  selection  rules  one  finds  that  the  deactivation  occurs  near  level 
v'  such  that 

v'.v'-l 

VNv_\  e'2weXehc/kT  (19) 

N0  ; 

1,2 

The  value  of  vNy  is  approximately  that  at  the  minimum  of  the  Treanor 
distribution  so  that  v'  may  be  obtained  in  terms  of  the  assumed  probability 
distributions.  Although  the  particle  flux  has  been  underestimated  in 
choosing  only  one  term  of  the  sum,  the  resultant  value  for  v'  is  quite 
accurate  because  of  the  rapid  variation  of  the  V-T  probability  with 
vibrational  level  as  shown  in  Fig.  1. 

The  quanta  source,  K,  necessary  to  maintain  9  is  found  to  be 
constant  for  the  different  V-T  probabilities  assumed  for  the  calculations. 
This  is  to  be  expected  since  there  is  no  direct  communication  between 
the  source  and  *Iie  deactivation  region,  and  the  only  requirement  on  the 
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source  is  that  it  maintain  the  temperature  in  the  lowest  levels.  Thus 

the  source  requirement  will  not  increase  until  V-T  probabilities  are 

sufficiently  large  to  affect  the  distribution  in  the  lowest  levels.  It  should 

be  emphasized  that  these  conclusions  regarding  the  limitation  to  the  quanta 

flux  and  the  constant  vNy  region  are  only  applicable  to  situations  in  which 

near  resonant  collisions  dominate.  In  fact,  it  is  shown  later  that  these 

calculations  are  invalid  if  the  effective  V-V  collisional  exchange  is  due 

to  nonre sonant  collisions  with  low  level  molecules. 

Figure  3  shows  vibrational  distribution  calculations  for  a  wide 

range  of  translational  temperatures  at  the  fixed  value  of  9  ^  =  2700°K. 

These  calculations  illustrate  the  effect  of  varying  the  relative  magnitude 

of  the  near  resonant  to  nonre  sonant  V-V  collisional  processes.  In  the 

low  level  region  controlled  by  detailed  balancing,  the  population  density 

of  the  higher  levels  decreases  with  increasing  translational  temperature 

as  seen  in  Eq.  (16).  The  typical  resonant  collisional  terms  in  the  quanta 

2  v,v+1 

flux  behave  like  N  P  ,  whilst  the  nonre  sonant  collisions  with  mole- 
v  v,v- 1  .  . 

v,v+l 

cules  in  the  lowest  levels  vary  as  NvN^  ^  .  Thus  for  high  vibrational 
levels,  the  collisions  with  the  more  densely  populated  lower  levels  become 
more  important  with  increasing  translational  temperature.  This  effect 
is  enhanced  by  the  increase  in  the  ratio  of  nonresonant  to  near  resonant 
probabilities  with  increasing  translational  temperature. 

The  results  for  the  three  lowest  translational  temperatures  all 
indicate  a  region  of  constant  vNv  implying  that  the  distribution  beyond  the 
Treanor  minimum  is  dominated  by  near  resonant  V-V  collisional  exchange. 
However,  at  the  higher  translational  temperatures  the  nonresonant 
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collisions  become  more  important  resulting  in  a  total  vibrational  inversion 
at  550°K.  Since  the  V-T  probability  increases  with  the  translational 
temperature,  collisional  deactivation  prevents  a  total  inversion  in  the 
distribution  at  700°K. 

The  possibility  of  a  total  inversion  can  be  readily  demonstrated 
when  nonresonant  collisions  dominate  i.  e.  when  the  most  effective  V-V 
collisions  are  those  with  the  highly  populated  low  vibrational  levels.  In 
this  case  the  leading  term  in  the  particle  flux  between  levels  v  and  v+1  is 

v,v+l  v+l,v 


F  ,  ,  ~  P 
v,v+l  ~ 


1.0 


N  N.  - 
v  1 


0,1 


Nv+1N0 


v+l.v 

and  this  is  equal  to  the  local  drain  which  is  approximately  Py  Nv+j* 
It  follows  that 


N, 


v+l,v 

V-T 


N 


v+1 


NjP 


v.v+l 


+  exp 


1.0 


/V1 

\ko; 


0  he 

2v  “exe  kT 


) 


kT  1  (20) 


For  vibrational  levels  beyond  the  Treanor  minimum,  the  above  expression 

v+l,v  v,v+l 

predicts  a  total  inversion  up  to  vibrational  levels  where  T^l^l  0 

is  O  (1).  The  approximations  in  the  above  derivation  are  invalid  near 
the  region  of  the  Treanor  minimum  because  of  errors  resulting  from 
small  differences. 

The  fact  that  the  present  calculations  for  CO  have  predicted  a 
small  total  inversion  for  only  a  very  limited  range  of  temperature  con¬ 
ditions  is  related  to  the  form  of  the  V-V  probability  distribution  used. 

v,v+l 

It  is  obvious  that  a  reduction  in  the  variation  of  P^  q  with  v  will  tend  to 
favor  the  nonresonant  collisions  and  should  lead  to  a  large  inversion. 
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This  is  verified  by  the  results  in  Fig.  4  where  we  have  used  the  artificial 

probability  distributions  as  shown  which  vary  the  ratio  of  near  resonant 

to  nonresonant  probabilities.  A  strong  inversion  is  evident  at  325° 

translational  temperature  for  the  flat  V-V  probability  distribution.  Even 

for  these  artificial  probability  distributions  near  resonant  collisions  are 

important  in  the  high  vibrational  levels.  Thus  the  source  is  not  linearly 

related  to  the  drain  and  the  inversion  is  limited  by  the  V-V  probability 

distribution.  It  should  be  emphasized  that  the  V-V  probabilities  shown 

in  Fig.  1  are  in  good  agreement  both  in  absolute  magnitude  and  v  depen- 

12 

dence  with  the  recent  experimental  data  of  Hancock  and  Smith.  Bray's 

2 

calculations  of  quasi  steady  state  distributions  in  Ng  also  predict  a 
total  inversion  after  the  minimum  of  the  Treanor  distribution  for  a  range 
of  values  of  0^  /T.  He  assumed  an  approximate  analytic  variation  of  the 
V-V  probabilities  with  vibrational  level  which  is  much  flatter  than  the 
present  probability  distribution  in  CO. 

B.  T-V  and  Radiative  Decay  in  Pure  CO 

Whereas  the  collisional  exchange  of  energy  between  the  vibrational 
and  translational  degrees  of  freedom  is  most  significant  at  the  higher 
vibrational  levels,  the  decay  of  vibrational  energy  due  to  radiation  is 
effective  in  all  vibrational  levels  for  an  optically  thin  gas.  In  the  present 
calculations  the  gas  is  assumed  to  be  optically  thin  for  all  vibrational- 
rotational  radiative  transitions.  This  is  valid  in  CO  for  pressures 
<  0.  1  Torr  at  room  temperature  over  a  characteristic  length  of  1  cm  for 
trapping  of  radiation.  These  conditions  are  typical  of  the  electrically 
pumped  CO  lasers.  The  results  of  these  calculations  are  relatively 
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VIBRATIONAL  POPULATION  (NORMALIZED) 


Fig.  4  CO  vibrational  distribution  functions,  neglecting  radiative  decay, 
for  0i  =  2700°K,  T  =  325°K  and  the  artificial  V-V  probability 
distributions  as  shown.  The  V-T  probability  distribution  used 
is  that  shown  in  Fig.  1. 
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insensitive  to  optical  opacity  in  the  first  few  vibrational  levels  because 
the  population  distribution  is  controlled  by  the  fixed  vibrational  tempera¬ 
ture  0  j' . 

At  low  pressure,  the  radiation  term  in  the  master  equation  (6) 
becomes  the  dominant  decay  term.  The  effect  of  the  radiative  vs.  V-T 
deactivation  is  shown  in  Fig.  5  as  a  function  of  the  CO  pressure.  For  the 
assumed  conditions  the  distribution  is  controlled  by  radiative  decay  at 
pressures  of  the  order  of  0.  1  Torr  or  less.  Because  of  the  distributed 
drain  in  energy  due  to  the  radiative  decay,  the  slope  of  the  vibrational 
distribution  no  longer  varies  linearly  with  v"  *  after  the  minimum  of  the 
Treanor  distribution. 

The  slope  of  the  distribution  in  the  high  levels  can  be  determined 
from  particle  flux  considerations.  At  low  pressures  where  the  radiative 
decay  is  the  dominant  deactivation  mechanism,  the  component  of  the  flux 
in  the  forward  direction  becomes  significantly  more  important  than  that 
in  the  backward  direction  and  it  follows  that 


N  v+l,v/  V*  v,v+1 


where  the  sum  may  be  evaluated  using  Eq.  (16)  for  the  first  few  levels. 

Nonresonant  collisions  are  more  effective  than  near  resonant  collisions 

because  of  the  low  population  densities  in  these  high  vibrational  levels. 

v.v+1 

Thus  the  slope  varies  approximately  inversely  with  Pj  q  .  It  is  evident 
from  the  above  result  that  the  high  level  population  distribution  depends 
critically  on  the  assumed  probability  distribution  when  the  radiation  decay 
dominates. 
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C.  Calculations  for  N2:  CO  Mixtures 

In  the  calculations  discussed  so  far  we  have  only  considered  pure 
CO  and  have  determined  the  distribution  for  a  fixed  temperature,  Oj  , 
in  the  lowest  levels.  Examples  of  the  more  general  case  in  which  the 
source  communicates  with  all  vibrational  levels  are  the  N^-CO  mixing 
experiments  in  which  the  CO  is  a  trace  species  and  the  "active"  acts 
as  an  infinite  source  of  vibrational  energy.  For  comparison  with  this  case 
we  have  assumed  a  fixed  vibrational  distribution  and  calculated  th*> 
steady  state  distribution  in  CO  as  a  result  of  energy  transfer  from  the 
to  the  CO  with  subsequent  intramode  energy  exchange  as  described  in 
the  preceding  sections. 

Because  of  the  large  energy  difference  between  the  first  vibrational 
levels  of  and  CO,  (188  cm  ^),  there  is  a  considerable  difference  in 
their  vibrational  temperatures.  For  example,  at  a  translational  temper¬ 
ature  of  32  5°K,the  CO  vibrational  temperature  of  2700°K  may  be  main¬ 
tained  by  a  vibrational  temperature  of  1  700°K  according  to  Eq.  (10). 

As  a  result  of  this  low  vibrational  temperature,  there  is  far  less 
intramode  pumping  of  vibrational  quanta  in  than  in  CO.  Thus,  the 
vibrational  distribution  is  quite  close  to  a  Boltzmann  distribution  at 
0  |  =  1700°K,  99%  of  the  molecules  and  92%  of  the  vibrational  quanta 
being  in  the  first  three  levels.  Therefore,  the  results  are  relatively 
insensitive  to  the  high  level  distribution.  For  the  present  calculations 
we  have  used  a  Treanor  distribution  up  to  its  minimum,  this  distribution 
being  necessary  to  maintain  V-V  equilibrium  between  the  anharmonic 
oscillators,  *  followed  by  a  constant  population  density  in  the  high  levels 
of  N2. 
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In  Fig.  6  a  comparison  is  made  between  the  vibrational  distribution 
calculated  for  0.  1  Torr  pure  CO  and  a  mixture  of  ^iCO  with  partial  pres¬ 
sures  of  2.  0  Torr  and  0.1  Torr,  respectively.  In  the  pure  CO  case  a 
vibrational  temperature  of  2  700°K  was  chosen  for  the  lowest  vibrational 
levc’s  whilst  the  mixture  calculations  assume  a  ^  distribution  which 
would  pump  CO  to  a  temperature  of  2  700°K  in  the  absence  of  any  energy 
drain.  In  the  mixture  calculations,  collisional  communication  is  allowed 
between  all  levels  of  and  CO.  The  mixture  calculation  leads  to  a 
vibrational  distribution  which  is  only  slightly  lower  than  the  distribution 
in  pure  CO  except  at  the  high  vibrational  levels  where  the  intramode 
pumping  cannot  keep  up  with  the  V-T  energy  loss. 

There  are  two  basic  reasons  for  the  difference  between  the  two 
distributions.  Firstly,  the  V-T  deactivation  in  the  high  levels  is  larger 
than  for  the  pure  CO  calculations  because  of  the  large  faction.  For 
the  present  calculations  we  have  assumed  that  the  probability  of  V-T 
deactivation  of  CO  by  ^  is  the  same  as  for  CO  by  CO  as  expected 
theoretically.  To  verify  the  effect  of  this  increased  V-T  drain,  we  have 
repeated  the  pure  CO  calculation  with  the  V-T  rate  increased  by  the  mole 
fraction  ratio.  The  result  is  shown  in  Fig.  6  and  the  effect  of  higher  V-T 
rates  is  only  evident  at  high  levels.  Secondly,  the  assumed  vibrational 
temperature  was  obtained  from  the  pumping  formula  (10)  without  taking 
the  energy  loss  mechanism  into  account.  It  is  this  small  but  finite  energy 
loss  which  reduces  the  population  distribution  and  the  vibrational  tempera¬ 
ture  0  j  at  the  lower  vibrational  levels.  Thus  to  achieve  a  particular 
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Fig.  6  CO  vibrational  distribution  functions  for  N2-CO  mixtures  with 
a  fixed  N2  distribution,  is  the  vibrational  temperature 
between  the  zeroth  and  first  levels  of  the  N2  distribution.  The 
partial  pressures  of  CO  and  N2  are  0.  1  and  2,  0  Torr, 
respectively,  and  T  =  325°K. 
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vibrational  temperature  in  the  lowest  levels  of  CO  requires  an  N2 
vibrational  temperature  somewhat  higher  than  predicted  by  Eq.  (10). 

For  a  harmonic  oscillator  model  one  can  readily  estimate  the 
correction  to  the  vibrational  temperature  as  a  function  of  the  ratio 
7y  when  radiation  is  the  dominant  loss  mechanism.  This 

correction  can  be  obtained  from  the  energy  relaxation  Eq.  (9).  For 
time  scales  short  compared  to  the  V-T  relaxation  time,  one  can  set 
the  last  two  terms  of  Eq.  (9)  equal  to  zero  in  the  steady  state  approxi¬ 
mation.  This  is  in  agreement  with  the  underlying  concept  in  this  work 
that  the  redistribution  of  vibrational  energy  occurs  in  time  scales  much 


shorter  than  the  collisional  deactivation,  i.  e.  7y  y«  Ty  j.  Then  one 
can  derive  the  following  relationship  for  the  vibrational  temperature 


in  terms  of  the  small  parameter  e 
vibrational  temperature  in  CO 


=  TV-v/^N2Trad 


and  the  required 


eCO  +  ®N2  “  0CO 

Tvco  tvn2  t 


In  1  +  e  (l-e‘0N2/TVN2) 


(22) 


For  the  sample  calculation  shown  in  Fig.  6  the  correction  to  the  N2 
vibrational  temperature  amounts  to  some  30°K.  The  corrected  CO 
distribution  is  shown  in  Fig.  6  and  is  v/ithin  a  few  percent,  of  the  pure 
CO  distribution  except  for  v  >  30  where  the  increased  V-T  deactivation 
is  important  as  described  above. 

It  is  evident  from  the  similarity  of  the  pure  CO  and  mixture 
calculations,  presented  in  Fig.  6,  that  the  vibrational  levels  of  N2 
communicate  predominantly  with  the  first  vibrational  level  of  CO  for 
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the  chosen  conditions  and  rates.  This  is  due  to  the  vibrational  energy 
defect  for  transitions  between  the  low  levels  of  Ng  and  higher  levels  in 
CO.  The  energy  defect  increases  with  vibrational  level  in  CO.  Con¬ 
versely  the  higher  levels  of  Ng  become  more  nearly  resonant  with  the 
v=l-*v=0  transition  of  CO  and  resonance  occurs  between  the  sixth  and 
seventh  levels  of  and  the  first  and  zeroth  level  of  CO.  Mixture  cal¬ 
culations  made  with  the  restriction  that  the  Ng  communicate  only  with 
the  first  and  zeroth  levels  of  CO  agreed  to  within  a  few  percent  with  the 
complete  communication  calculation  for  the  assumed  V-V  rates.  It 
should  be  noted  that  a  sufficiently  large  increase  in  the  I^-CO  V-V  rate 
would  result  in  an  increase  in  the  high  level  population  densities  in  CO 
as  a  result  of  some  direct  communication  with  the  low  levels  in  N£. 
However,  a  decrease  in  V-V  rate  with  suitable  correction  to  the  Ng 
temperature  by  means  of  Eq.  (22)  would  not  change  the  CO  distribution 
since  the  would  then  only  communicate  with  the  lowest  levels  of  CO. 

Thus,  the  mixture  calculations  show  that  for  typical  ^-CO 
mixtures  the  distribution  in  CO  is  characterized  mainly  by  the  low  level 
temperature  0^  .  The  nonresonance  in  the  V-V  energy  exchange  between 
and  CO  results  in  the  CO  intramode  exchange  being  much  faster  than 
the  effective  intermode  exchange  when  corrected  by  the  mole  fraction 
ratio. 

Comparison  with  Experimental  Data 

The  currently  available  experimental  data  for  the  vibrational 
distribution  in  CO  have  been  obtained  in  electrical  N^-CO  mixing 
experiments.  The  Ng  is  first  made  vibrationally  active  by  passage 
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through  a  RF  discharge  and  subsequently  undergoes  V-V  collisional 

exchange  with  CO  which  is  mixed  in  small  concentrations  into  the  N_ 

15  17 

stream.  In  the  experiments  by  Joeckle  and  by  Horn  and  Oettinger, 

the  vibrational  distribution  in  the  CO  is  derived  from  the  overtone 

1 8 

emission  spectrum.  Legay  obtained  data  for  the  vibrational  distribu¬ 
tion  from  v=7  to  v=17  from  gain  measurements  in  an  electrical  ^-CO 
mixing  laser.  In  order  to  compare  the  experimental  results  with  the 
steady  state  predictions,  it  must  be  assumed  that  quasi  V-V  equilibrium 
has  been  reached  in  the  experiment.  This  was  experimentally  verified 
by  Horn  and  Oettinger  and  is  almost  certainly  true  in  the  other  two 
experiments  in  which  the  time  after  mixing  is  long  compared  with  the 
CO  V-V  relaxation  time. 

All  the  calculations  for  comparison  v/ith  the  experimental 

vibrational  distribution  data  have  been  obtained  for  pure  CO  with  the 

>!< 

appropriate  8j  and  due  allowance  made  for  the  increased  V-T  deacti¬ 
vation  due  to  the  N£.  The  pure  CO  calculation  is  justified  in  view  of 
the  mixture  calculations  in  the  preceding  section  which  were  made  for 
similar  mixtures  of  and  CO.  Furthermore,  it  is  assumed  chat 
metastable  electronic  states  in  the  excited  play  no  role  in  the 
excitation  of  the  CO.  Figure  7  shows  a  comparison  of  the  experimental 
data  of  Joeckle  and  the  numerical  results.  Since  8^  is  not  measured 

>jc 

in  the  experiment,  we  have  chosen  8^  to  provide  the  best  fit  of  a  Treanor 
distribution  with  the  low  level  data  making  use  of  the  experimentally 
determined  translational  temperature.  Two  numerical  results  are 
shown,  the  second  being  derived  with  twice  the  nominal  V-V  probability 
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C  6060 


VIBRATIONAL  LEVEL 


.  7  Comparison  of  experimental  CO  distribution  (Ref.  15)  with 

calculations  for  0i  k  =  2500°K,  T  =  350°K  and  Pqq  =  ^8  Torr 

The  upper  theoretical  curve  illustrates  the  effect  of  increasing 
the  V-V  probabilities  by  a  factor  of  2. 


to  demonstrate  the  sensitivity  of  low  pressure  data  to  the  V-V  rates. 

It  is  evident  from  the  master  equation  that  a  change  in  the  magnitude  of 
the  V-V  probabilities  is  equivalent  to  a  change  in  the  CO  pressure  by 
the  same  amount.  This  can  be  significant  at  low  pressures  when  the 
radiative  decay  is  the  dominant  loss  mechanism  as  in  the  present  results. 

The  data  of  Horn  and  Oettinger  are  shown  in  Fig.  8  and  correspond 
to  lower  pressures  in  CO  and  much  higher  vibrational  temperatures. 

jjc 

The  numerical  calculations  have  been  made  for  two  values  of  9^  for 
each  pressure,  the  lower  value  corresponding  to  that  determined  by 
Ho^  i  and  Oettinger  in  curve  fitting  a  Treanor  distribution  to  the  first 
four  vibrational  levels.  Their  values  of  9^  are  a  lower  estimate  because 
of  the  low  level  deactivation  due  to  the  strong  effect  of  radiative  decay 

5|C 

at  the  low  pressures.  This  underestimate  of  9^  is  evident  in  the 

❖ 

improved  fit  of  the  data  for  the  calculations  at  the  higher  values  of  9^  . 
These  higher  values  of  9j  were  chosen  to  demonstrate  the  effect  of 
low  level  radiative  decay  and  are  not  necessarily  the  optimum  values. 

The  sensitivity  of  the  calculations  to  the  V-V  probabilities  iE  somewhat 

$ 

less  than  in  the  comparison  in  Fig.  7  because  of  the  higher  values  of  9j  . 

Jje  q 

The  experiment  at  9j  =3500  K  has  the  largest  mole  fraction  ratio  and 
the  highest  temperature,  approximately  2060°K  in  the  lowest  levels, 
of  all  the  data  considered.  A  mixture  calculation  for  this  extreme  case 
agreed  with  the  pure  CO  distribution  calculation  in  the  lower  levels  and 
was  less  than  a  factor  of  two  larger  in  the  higher  levels. 

Finally,  we  have  the  comparison  with  the  measurements  by  Legay 
presented  in  Fig.  9.  These  measurements  were  made  at  a  much  higher 
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VIBRATIONAL  POPULATION  (NORMALIZED) 


Fig.  8  Comparison  of  experimental  CO  distribution  (Ref.  17)  with 

calculations  for  T  =  350°K.  The  lower  curves  are  for  values 
of  deduced  in  Ref.  17,  and  the  upper  curves  show  the 
effect  of  a  small  variation  in  01*. 
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VIBRATIONAL  POPULATION  (NORMALIZED) 


O  EXPERIMENT 
—  THEORY 


VIBRATIONAL  LEVEL 


CO  distribution  calculation  compared  with  the  experimental 
data  of  Ref.  18  for  ©i*  =  3000°K,  T  =  325°K  and  PCO  =0-3  Torr. 
The  experimental  data  are  limited  to  v>7  as  they  were  deduced 
from  gain  measurements  in  a  N2-CO  laser  which  only  operated 
above  v=6. 


C  O  pressure  than  the  previous  data  and  the  calculations  do  not  predict 
the  high  level  decay  which  was  observed  experimentally.  Although  there 
is  a  possible  large  error  in  the  estimated  value  of  0  j  ,  this  will  not 
greatly  affect  the  high  level  population  density  variation  because  radiative 
effects  are  expected  to  be  relatively  small  at  0.  3  Torr  as  indicated  in 
Fig.  5.  There  is  no  obvious  explanation  for  the  discrepancy  between 
theory  and  experiment  at  the  high  vibrational  levels.  One  possible 
reason  for  the  apparent  discrepancy  in  the  high  level  data  is  the  deactiv¬ 
ation  of  these  levels  by  O^  v/hich  is  present  in  the  laser  to  maintain  an 
equilibrium  concentration  in  CO.  The  high  vibrational  levels  in  CO  are 
resonant  with  the  v=  1  to  v=0  transition  in  O^.  The  experimental  technique 
may  also  result  in  small  inconsistencies  in  the  data  because  it  requires 
separate  gain  measurements  for  each  pair  of  vibrational  levels  unlike 
the  overtone  spectrum  which  encompasses  all  vibrational  levels.  A 
reduction  in  the  V-V  rates  would  reduce  the  high  level  population,  but 
this  reduction  is  not  suggested  by  the  comparison  with  the  other  experi¬ 
mental  data.  It  is  also  possible  that  the  high  level  V-T  probabilities 
have  been  underestimated  in  the  present  calculations. 

We  believe  that  the  comparisons  shown  in  Figs.  7  to  9  indicate 
reasonable  agreement  of  the  experimental  data  with  the  analytic  model 
of  vibrational  energy  transfer  in  CO  and  I^-CO  mixtures. 
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IV.  CONCLUSION 


A  model  has  been  developed  to  calculate  the  steady  state  vibrational 
distribution  of  a  diatomic  gas  under  conditions  of  thermal  nonequilibrium. 
Calculations  have  been  made  for  the  case  of  a  binary  gas  mixture  with 
one  gas  acting  as  an  infinite  source  of  energy  for  the  other.  It  is  shown 
that  under  suitable  conditions  the  distributions  may  be  well  approximated 
by  the  assumption  that  the  source  of  energy  which  maintains  the  distri¬ 
bution  communicates  only  with  the  zeroth  and  first  vibrational  level. 

The  model  can  readily  be  extended  to  predict  the  distributions  in  a  binary 

1  9 

mixture  excited  by  an  external  source.  An  example  of  interest  is  the 
directly  excited  N^-CO  laser.  The  calculations  include  both  collisional 
deactivation  and  radiative  decay  to  account  for  the  energy  loss  from  the 
system.  These  decay  processes  prevent  a  large  inversion  in  the  upper 
vibrational  levels  which  might  otherwise  occur  as  a  result  of  intramode 
vibrational  energy  pumping  due  to  the  molecular  anharmonicity.  It  is 
shown  that  radiative  decay  may  be  the  dominant  loss  mechanism  in  low 
pressure  N^-CO  experiments. 

The  numerical  calculations  are  found  to  be  in  reasonably  good 
agreement  with  the  available  experimental  data  all  of  which  have  been 
obtained  in  electrically  excited  ^-CO  mixing  sys'ems.  It  is  evident 
from  the  comparison  with  the  experimental  data  that  more  detailed 
kinetic  rate  information  is  needed  particularly  for  the  high  level  colli¬ 
sional  deactivation  processes.  Within  that  limitation  the  present  model 
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gives  a  satisfactory  account  of  the  experimentally  observed  non  Boltzmann 
distributions  based  solely  on  collisional  exchange  of  vibrational  energy 
and  radiative  decay. 

On  the  basis  of  existing  theories  for  the  probability  of  collisional 
exchange,  it  would  appear  that  near  resonant  intramode  pumping  will  be 
an  important  if  not  the  dominant  mechanism  for  energy  pumping.  When 
near  resonant  collisions  are  the  dominant  intramode  pumping  mechanism, 
it  has  been  shown  that  a  total  inversion  is  not  possible  unless  the  source 
can  communicate  directly  with  high  vibrational  levels  by  intermode 
pumping. 

Under  conditions  where  the  near  resonant  collisions  dominate 
and  radiative  decay  is  not  important,  some  simple  criteria  have  been 
derived  to  predict  the  distribution  function.  These  predictions  are  based 
on  knowledge  of  the  translational  temperature  and  the  low  level  vibra¬ 
tional  temperature  which  are  more  readily  accessible  to  experimental 
measurement  than  the  high  level  distributions.  Similar  prediction 
criteria  have  been  derived  for  the  case  when  radiative  decay  is  the 
dominant  energy  loss. 
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